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Abstract 

The presence of extra scalar singlets is a feature of several motivated extensions of the 
Standard Model, and the mixing of such a singlet with the Higgs boson is allowed to be 
quite large by current experiments. In this paper we perform a thorough phenomenological 
study of this possibility. We consider both direct and indirect searches, and we quantify the 
current constraints as well as the prospects for future hadron and lepton machines - from 
the forthcoming LHC run up to a futuristic 100 TeV proton-proton collider. The direct 
reaches are obtained extrapolating the current limits with a technique that we discuss and 
check with various tests. We find a strong complementarity between direct and indirect 
searches, with the former dominating for lower values of the singlet mass. We also find that 
the trilinear Higgs coupling can have sizeable deviations from its Standard Model value, a 
fact for which we provide an analytical understanding. The results are first presented in a 
general scalar singlet extension of the Standard Model, taking advantage of the very small 
number of parameters relevant for the phenomenology. Finally, we specify the same analysis 
to a few most natural models, i.e. the Next-to-Minimal Supersymmetric Standard Model, 
Twin Higgs and Composite Higgs. 
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1 Introduction 


Is the Higgs boson 0@ the only light scalar particle, or is it rather the low-hanging fruit of 
a richer scalar sector of Nature? This question is of paramount importance for the field of 
particle physics. Discovering another scalar close-by would in fact have deep implications for 
our understanding of the concept of naturalness 3 , challenging frameworks that do not aim at 
removing the apparent tuning of the weak scale 4 . 

Among the many conceivable scalar extensions of the Standard Model (SM), arguably the 
simplest one consists in the addition of a singlet under the SM gauge group which mixes with 
the already discovered Higgs. This possibility is ubiquitous in many new physics scenarios, most 
notably in two of the most successful proposals to address the hierarchy problem of the Fermi 
scale: the Next-to-Minimal Supersymmetric Standard Model (NMSSM, see for a review), 
and Twin Higgs [^. It is also interesting that, in the context of “finite naturalness” [^, an 
extra singlet-like scalar would also allow to dynamically generate the weak and Dark Matter 
scales 1^. Moreover, and independently of the hierarchy problem, the same addition to the SM 
among the simples possibilities to achieve electroweak baryogengesis [9-11 


IS 


From an experimental point of view, the 8 TeV LHC run has given us a knowledge of 
the main Higgs properties to a 20% accuracy [12, 13 , and no direct evidence for new physics 


whatsoever. Despite that, scenarios exist that are poorly constrained by present data, an extra 
scalar singlet mixed with the Higgs being a paradigmatic example of this case. 

The above considerations strongly motivate a thorough experimental quest for such a parti¬ 
cle. The goal of this paper is to provide a general and solid strategy to perform this hunt both 
at current and future machines. 

Several studies in this direction have been performed in the literature (see e.g. the recent 
]14| - [T7| ), however they have mostly restricted their attention to some particular direct and/or 
indirect search, to some particular model, or to some particular machine. In this paper our aim is 
to be more comprehensive, building a picture that allows to grasp the general phenomenological 
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features, yet in an exhaustive way. In order to do so, we develop an analytical understanding 
of what the parameters most relevant for the phenomenology are, freeing us from the necessity 
of numerical scans. In this way we are able to identify the most promising ways to probe this 
scenario, thus providing key information for current experiments, as well as for planning future 
ones. 

In particular, we will consider several experimental programs besides the second run of the 
LHC that has just started. Our study will include the high luminosity phase of the LHC itself 
(HL-LHC [^), a higher energy phase (HE-LHC), linear (ILC [19| and CLIC [20[|21] ) as well 
as circular (FCC-ee and CEPC [^) electron-positron colliders, up to a futuristic 100 TeV 
proton-proton collider (ECC-hh |24|). 

The natural approach to study such an extended Higgs sector relies on two complementary 
strategies: 

o Indirect searches: modification of the couplings of the 125 GeV Higgs boson; 

o Direct searches: production and detection of the heavier scalar. 


If on one hand the impact of an improvement in precision on the Higgs signal strengths is quite 
easy to estimate, on the other hand the study of the reach of direct searches requires generally 
an amount of numerical analysis. In this paper, however, we rely on a technique introduced 


m 


25 and developed in 26 to extrapolate the reach of future colliders from the existing limits 


of the first LHC run. This will help to get an intuition of the capability of the new experiments, 
avoiding the use of detailed collider simulations, and it will allow for an easy comparison between 
the relevance of precision physics and direct searches. 

The paper is organised as follows. In sectionwe introduce the model and present the basic 
phenomenological consequences, showing also what are the simplifications that can be made. In 
section we review the extrapolation procedure presented in 25,26 discussing its virtues and 
limitations, and we present our projections for the sensitivity of direct searches in the relevant 
channels at future machines. These results have a general validity that goes beyond the context 
of the singlet scenario analysed in this paper. We then apply the bounds to the case of generic 
singlet scalar coupled to the SM. Sections and are dedicated to the NMSSM, and Twin 
and Composite Higgs, respectively, and there we project present and future bounds into the 
relevant parameter space of the two models. We conclude in section Technical details about 
the extrapolation procedure are given in appendix 


2 Phenomenology of a scalar singlet 

2.1 A convenient description 

We add to the SM particle content a scalar field S, singlet under the SM gauge group. The 
quantum numbers of the real component of S, on which we will focus, allow it to mix with the 
SM Higgs boson, and we dub the mixing angle 7 . Let h and (p be the two neutral CP-even 


propagating degrees of freedom, h being the one observed at the LHC, and ruh = 125.1 GeV 27 
and their physical masses. In terms of the gauge eigenstates 


vr' 


H = I V + + Z7r° 

V2 


0 


5 = 


Vs + s' 


the mass eigenstates read 

h = hP cos 7 -|- s° sin 7 , (p = —hP sin 7 -|- s*^ cos 7 , 

where Vg is the vacuum expectation value (vev) of the singlet field, and v ~ 246 GeV. 


( 2 . 1 ) 


( 2 . 2 ) 
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All the h and cj) couplings to SM fermions and vectors are fixed, dne to the singlet nature 
of S, only by the mixing angle 7 . They can be written in terms of the couplings of a standard 
Higgs boson of the same mass as 


9hff 


9hVV 

9WV 


= c 


'75 


9<t>ff ^ 9<j>VV 

„SM „SM 

9hff 9hVV 


— S. 


1 > 


(2.3) 


where we introduce the notation sq = sin0, cg = cos 9. Notice that the branching ratios of h 
are not modified with respect to their SM values, since all the couplings are rescaled by the 
same factor; the only observable deviation from the SM is therefore the reduced production 
cross-section. Concerning cj), its production cross-section is simply the one of a standard 

Higgs boson of mass rescaled by Its branching ratios are again those of a SM Higgs 
boson of the same mass below the kinematic threshold = 2 m/i, where the width 
can in principle become sizeable. All these considerations can be modified by the presence of 
operators of dimension greater than four in a strongly coupled theory, but in the following we 
will consider only weakly coupled scenarios. 

The phenomenology of the two scalars can then be summarised as 

fih = X fisM, (2.4) 

9(l>^vvjf = s'^ X usMirritj,) X (1 - BR^p^hh ), (2.5) 

X (Tg]vi(?^(/>) X , (2-6) 


where are respectively the signal strengths of h and (f) in any SM channel (unless specified), 
^SM is the signal strength of a SM Higgs boson, and cjsm is its production cross-section. The 
measurement of the h signal strengths, and hence of its couplings, puts a model independent 
constraint on the mixing angle 7 . This in turn limits the magnitude of possible direct signals of 
(/), which for large enough rriyp decays mainly into a pair of vector bosons, W^W~ and ZZ, and 
into hh if the relative branching ratio is large enough. The only relevant fermionic decay mode 
of (/> is a pair of top quarks, which is however subdominant with respect to the previous ones. 

One reads in expressions (2.5) and (2.6) that the decays of both scalars are completely 
determined by just three parameters: 


o the mass of the singlet-like scalar m^, 
o the mixing angle 7 , 


o the branching ratio BRpy^hh- 

We now discuss them in order to understand better and simplify the description of the phe¬ 
nomenology, starting with BRyp^hh. 

For nifp ^ rnw, one knows that decays of cj) into vector bosons are determined by their 
Goldstone nature, and that they dominate over all the other SM decay modes. This implies the 
asymptotic relation 

BRpy^hh = BR<^_s.zz = - BRpy^ww = - , > mw- (2.7) 

For large enough the number of free parameters drops then from three to two, further 
simplifying the description of the model. It is natural to choose as one of them. 

Concerning the other free parameter, notice that for any given scalar potential, the mixing 
angle 7 between the two CP-even scalars can be expressed in the form 


• 2 ^hh ~ '^h 

sm^7 = —-ft. 

' 


( 2 . 8 ) 
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where is the ( 1 , 1 ) element of the 2 x 2 mass matrix of the scalar system before 
diagonalisation, in the {h^, s^) basis. Notice that one has three possibilities to obtain a relation 
between 7 , the physical masses and the parameters of the Lagrangian: one possibility for each 
independent entry of the mass matrix M^. We find it convenient to choose because it 

is always of the order of the electroweak scale v, despite its explicit form as a function of the 
Lagrangian parameters depends on the model under consideration. While this statement will 
be made more explicit in the discussion of the potential in section | 2 . 2 [ it can be understood 
qualitatively by noting that, in general, it applies to all cases where the only field responsible 
for EW symmetry breaking (EWSB) is the one associated with that single entry of This 

property holds in particular for the doublet-singlet system under consideration, since the scalar 
field S does not play any role in EWSB. We conclude that the use of as a free parameter 
to describe the phenomenology, instead of sin^ 7 , is a convenient choice. 


2.2 Potential and couplings among the Higgses 

The discussion performed so far allows to treat the decay (/>—)• /i/i only in an approximate 
manner. Moreover, it has nothing to say about possible deviations of the trilinear Higgs coupling 
9hhh from its SM value. It is known that such deviations are allowed to be quite large, at least 
in some specific models containing singlets [28] . It would then be interesting to compare them 
with the other indirect signals, possibly in a most general way. 

With this in mind, we write down the most general form of a renormalisable scalar potential 
involving a doublet H and a singlet S as 

V{H, S) = f4j\H\‘^ + Xh \H\^ + Xhs S^\H\^ + an S\H\^ + l^s + A 5 (2.9) 


By deriving the potential with respect to S, one sees that the condition = 0 for the singlet 


vev requires an = 0. This further implies sin 7 = 0, as one can infer from (2.9). This is not the 


case we are interested in, and hence we assume Us 7 ^ 0 in the rest of our discussion. Since 7 ^ 0 
spontaneously breaks the symmetry S —)■ —S', we assume in general that also an and as are 
different from zero. A non zero positive A 5 is then a necessary condition for the potential to be 
bounded from below (neglecting runnings at higher energies that could change this limitation, 
see e.g. (Sl). 


We can then substitute a subset of the seven parameters of the potential in (2.9), with 


the vevs v and Vg and with the physical masses nih and m^. We start by removing the mass 
parameters fin and /X 5 , in favour of the vevs. In this way one realises that the following relation 
holds 


Mlh = 2A//u^ 


( 2 . 10 ) 


substantiating the discussion of the previous section. We then proceed to exchange the last two 
dimensionful parameters of (2.9), an and as, in favour of the physical masses nih and The 


most general addition of a scalar singlet to the SM can then be described in terms of five free 
parameters: 

Mhh, Vs, Xhs, ■ (2-11) 

The mixing angle 7 , which determines most of the phenomenology of the model, is a function 


of only two of them, see (2.8). The couplings among the scalars, on the other hand, explicitly 


depend also on the other three parameters. They read 


^The only dimensionful parameter of the SM potential is a trivial example of the above property: its size is 
related to v by an order one factor, the Higgs quartic coupling. 
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(2.13) 


where 7 is as in (2.8) and is the SM value of the triple Higgs coupling. One would like 
to know how Vs, Xhs and Xs impact the phenomenology, compared to m,ji and Mhh'- if all the 
parameters of ( 2 . 11 ) were equally important, as one could expect in principle from ( 2 . 12 ) and 


(2.13), the description of this part of the phenomenology would be somehow obscured. We show 


here that this is not the case, and that the impact of Vs actually dominates over the one of Xhs 
and Xs in most of the relevant parameter space. We do so by performing a large m^ expansion 
of the interesting phenomenological quantities, and ghhh/gf^h-, which reads 


BR 


(p^hh 


1 3 V 

A A Vs m^ 


, 3 - ml 

8 U 2 ml 


2m\r + m| + 2ml - - 2 Xhsv^ 


rrij 


1 3 V 
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9hhh 
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ml 
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1 ML - m 
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I I 3ml - 2 X hsv 

2 ml 



(2.15) 


Here we are assuming that and Vs are of order of the weak scale, and we collectively denote 
by 0{v‘^ jml) orders of Xhsv‘^/ ml, Xsvl/ml, and We do not keep track of the explicit 

dependence on the couplings in the higher order terms, since we are considering only weakly 
coupled scenarios where these couplings are never much greater than one. A few remarks are 
in order: 


o both quantities have the correct asymptotic behaviour if the singlet is decoupled; 

o at first order in v/m^, only Vs determines and ghhh/gLth (besides and Mhh), 

and at second order only Xhs enters the expressions. We have thus obtained an analytical 
understanding of which extra parameters will impact more on the phenomenology; 


o if v/vs = 0{m^/v) or bigger, then both expansions (2.14) and (2.15) are inconsistent. 


In the following, we will always use the complete expressions (2.12), (2.13) in our phenomenolog¬ 
ical studies. We notice that, choosing a very small value for Vs/v, the width of the singlet-like 
scalar becomes larger than its mass, and the triple Higgs coupling grows as jv/, thus making 
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Qhhh 

9hfh 


1(7 reach in 


LHC8 

LHC14 

HL-LHC 

HE-LHC 

FCC-hh 

ILC 

ILC-up 

CLIC 

CEPC 

FCC-ee 


0.2 

0.08-0.12 

4-8x10-2 


2 x 10-2 

4x10-3 

2-3x10-3 

2x10-3 

1x10-3 


0.5 

0.2 

0.08 

0.21-0.83 

0.13-0.46 

0 . 1 - 0.21 


Table 1. Current and indicative expected precisions on s2, and on the trilinear Higgs coupling ghhh- 


The values for s2 correspond to the most precise coupling (always ZZ and/or WW). All the numbers 


are taken from 32 (see 23 for CEPC). 


the model effectively strongly coupled. Furthermore, as also grows for small introducing a 
hierarchy among the dimensionful parameters in the potential. 

We conclude this section observing that the minimum {{\H\), (S)) = (v,Vs)/V^ is always 
the lowest one, for a given choice of the parameters. We verified this in a range that widely 
covers our phenomenological study. For a more detailed study of the theoretical constraints on 
a singlet-Higgs potential, see (^ . 

2.3 Indirect signals: triple Higgs conpling vs. signal strengths 

We have now developed enough tools to study the impact of indirect searches in the scenario 
under consideration. Before exposing quantitative predictions, it is important to know the 
present and future experimental situation. We thus begin this section with an overview of 
the current exclusions, as well as the expectations for future sensitivities in precision Higgs 
measurements. 


Current and future sensitivities to Higgs couplings. A global fit to all the Higgs signal 
strengths as measured at the 8 TeV LHC constrains sin^y < 0.23 at 95% C.L. 30,31 . Con¬ 


cerning future projections, we rely on the comprehensive summary performed at Snowmass in 
2013 32 , and we refer the interested reader to it for more detailsj^ 


We report in table the expected Icr precisions in coming from the measurement of a 
single coupling; for comparison we provide also the current best precision in a single channel 
from the LHC8 data. In the case of the mixing with a singlet, these numbers provide at least 
a good indication of the possible outcome of a global fit. Notice that this is not the case for 
more complex Higgs sectors, like when the Higgs mixes with another doublet (e.g. MSSM), 
where a global fit is usually much more constraining (see e.g. [28]). The hadronic collider 
configurations that we consider are LHC at 14 TeV with 300 fb^ (LHC14) and 3000 fb-^ 
(HL-LHC) of integrated luminosity. Assuming one will measure SM central values of the Higgs 
signal strengths, the expected precision will be in the ballpark of 10% and 5% respectively. The 
ranges reported in table come from different assumptions on the future theory errors [33| , and 
from different systematics scenarios |34| . Future projections for higher energy hadronic colliders 
are currently unavailable, needless to say they would constitute an essential ingredient in the 
present analysis. For leptonic machines, we report the expected precisions for the ILC and its 
luminosity upgrade (ILC-up), for CLIC and for circular colliders such as FCC-ee or CEPC. For 


More recent results, when available (see e.g. 13 33 ), substantially agree with the numbers given here. 
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Figure 1. Exclusion at 95% C.L. from current Higgs couplings measurements at the LHC8 (coloured 
region), and deviation in the Higgs signal strengths (coloured lines). In both plots the ratio of the 
trilinear Higgs coupling to its SM value is drawn for Vg = 250 GeV (dashed black). Grey: unphysical 
parameters. 


CLIC we report two numbers, the first assuming a data taking up to 1.4 TeV with 1.5 fb“^, the 
second adding 2 fb“^ at 3 TeV. Assuming SM central values, could be constrained from the 
per-cent down to the per-mille level, FCC-ee being the most sensitive collider in this respect. 
Concerning the trilinear Higgs coupling, for definiteness we stick again to the review 32 , and 


refer the interested reader to it for details. The expected Icr precisions achievable at future 
machines are also reported in table The measurement of the trilinear Higgs coupling benefits 
from higher energy collisions, thus the projection for FCC-ee is not reported. The benchmarks 
we consider for hadronic (pp) colliders are a 33 TeV one (HE-LHC) and a 100 TeV one (FCC- 
hh), both with 3 ab“^ of integrated luminosity. The two numbers given for ILC and ILC-up 
correspond to the ILC500 and ILCIOOO cases, those given for CLIC have the same meaning 
as for the signal strengths, and assume electron beam polarisation (for unpolarised beams, the 
worst expected precision would increase to 0.28). Again we warn the reader that such numbers 
are to be taken with a grain of salt, possibly even more than in the case of signal strengths 
In any case it looks very challenging to push the sensitivity in ghhh below the ~10% level. 


32 


Results and interplay of different measurements. The predictions for and Qhhh are 
shown in figures [T] and The shaded pink region is excluded by a combined fit to Higgs coupling 
measurements 30,31 , while continuous pink lines display some values of s^. Given that the 


estimates for future sensitivities are only indicative, and that those for some machines (like 
FCC-hh) are not yet available, we prefer not to present lines corresponding to specific future 
colliders. This will make it easier to visualise the reach due both to new precision assessments, 
and to new data. With this in mind and to avoid clutter, in the left-hand plots we focus on 
a region of parameters interesting for LHC and its upgrades, and in the right-hand ones we 
zoom into an area more interesting for future leptonic machines. One can see that reaching a 
per-mille precision on the Higgs signal strengths would constitute an impressive probe of this 


scenario. This will become even clearer when compared with direct searches in section 3.3 
We show the predictions for ghhh/dhhk dashed black lines. As discussed in the previous 
section, this quantity is not univocally fixed in terms of and M^hi so we show it for the 
benchmark values Vg = —75 GeV and Vg = 250 GeV. One can see that the allowed parameter 
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Figure 2. Same as figurebut with Vs = —75 GeV. 

space can accommodate very large deviations from one in ghkh/gf^h^ with both signs and up 
to an order of magnitude. From (2.15) one infers that the large size of such deviations is quite 
general, in other words it is not typical of the particular values chosen for Vg, at least until v/vg 
does not become too small. We explicitly verihed this by exploring other values of Vg not too 
far from v. We further check numerically that variations of both Xhs and Xs have a minor 
impact in the size of ghhh/gf^h^ confirming the insight from our analytical formula (2.15) (the 
plots being for Xhs = = !)• Finally, another check we perform is that, in the whole interval 

of parameters that we consider, one never has F,^ > rriij,. This would be the case e.g. if we have 
chosen a value of Vg closer to zero. 

An important message of figures [l] and and table [T] is that, at each stage of the future 
experimental program (except for FCC-ee and CEPC), the trilinear Higgs coupling and the 
Higgs signal strengths will be both crucially important as probes of the model. Moreover, 
under some circumstances, ghhh could even be the only observable where deviations will show 
upQ Perhaps more important than that, in the currently allowed parameter space one could 
have deviations which are so large, to induce an increase of an order of magnitude (or more) of 
double Higgs production (see e.g. [^), likely already in the reach of current machines. 

The presence of so large deviations in the trilinear Higgs coupling is a rather peculiar prop¬ 
erty: the opposite is in fact true in many other models (see e.g. [^), for example when the 
Higgs mixes with another doublet. Our study shows an explicit and rather general model where 
such large deviations are actually typical. It then motivates a further effort in the community 
to measure the Higgs self-coupling, already with the available data, and to assess the relative 
sensitivity of future machines. 


Electroweak precision tests. In the model under consideration the only contribution to the 
electro-weak precision tests (EWPT) is due to the IR-logarithms in the S and T parameters, 
that can be estimated as fSSl 


S 


a 


ml 

log-4, 


mf 


f = - 


3a 

16 TT cl, 



n 


(2.16) 


where a is the hne-structure constant and and are the sine and cosine of the Weinberg 
angle. S and T are proportional to the reduction of the h couplings to the weak bosons, i.e. to 


^This can be expected in a bottom-up perspective from the EFT analysis of 
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the mixing angle. They also vanish when being the mixing unphysical in this limit. 

EWPT do not have a strong impact on the picture under consideration, unless does not go 
to zero when is decoupled. For this to be possible, one needs M^h to be proportional to 


rriif) (see (2.8)): while this is not the case in Supersymmetry, Twin and Composite Higgs models 
realise such a behaviour. We will then return on the EWPT in section! 


3 Direct searches 

Searches for scalar resonances decaying into a pair of SM particles have been performed by the 
ATLAS and CMS experiments in a large variety of channels. Being based on the full dataset 
of the first LHC run, many of these searches are becoming sensitive to cross-sections that are 
significantly smaller than the ones of a standard Higgs boson. Therefore they already set relevant 
constraints for scalar singlets. With the increasing center-of-mass energy and luminosity, all 
these bounds will strengthen further. In this section we discuss the searches that are relevant 
for our purposes, and present extrapolations of their limits to future colliders. 

3.1 Extrapolation of the bounds 

In order to get a projection of the mass and cross-section reach of future collider experiments, we 
extrapolate the expected limits from the 8 TeV LHC relying on a parton luminosity rescaling, 
as introduced in and developed in 26 . This method, which we are going to recall here, 


is subject to a number of assumptions that are potentially not satisfied in a generic search. It 
nevertheless constitutes a quick and efficient tool to get a reasonable estimate of the masses and 
cross-sections that could be probed, avoiding to perform a full Monte-Carlo simulation of the 
experimental analysis. 

The main underlying idea is that the number of signal events that can be excluded with a 
certain confidence level at a given collider, for a given mass of the resonance, depends exclusively 
on the corresponding number of background events Nb- In other words, starting from the LHC 
bounds, it is possible to obtain a projection for the exclusions at a different collider, as a 
function of the resonance mass, studying the scaling of the background with the collider energy 
and luminosity. 


Extrapolation procedure. 

mass, mo 


Concretely, one proceeds as follows 26 


Take a value of the 

and a corresponding value of the cross-section times branching ratio, [cr x BR]o, 


that can be excluded at 95% C.L. at a collider with center-of-mass energy yTo and integrated 
luminosity Lq. It is possible to find a new value of the mass, m, for which the same number 
of background events is produced at a different collider with energy ^/s and luminosity L, by 
solving 

NB{m;s,L) = NB{mo;so,Lo). (3.1) 

The number of signal events that can be excluded for the new mass m with the new collider 
will then simply be equal to the one of the original search, Ns{m] s, L) = Ns{mo] so,Lq). 

If the signal acceptances and efficiencies do not change significantly between the two ex¬ 
periments - an assumption which we shall discuss further below - one can finally get the new 
value of the excluded cross-section dividing the number of signal events by the luminosity L, 


and therefore, imposing (3.1), 


[a X BR](m;s,L) = ^[a x BR]o(mo; sq, To)- 


(3.2) 


The extrapolation of the backgrounds to higher energies and luminosities, which is needed to 


determine the mass m from (3.1), is performed in the following way. The number of background 
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events of a certain process can be written as 


NB{m; s, L) (X L ■ 

{i,j}' 


, ^ d/ijj (s, s) ^ 
ds-- CTijis), 


(3.3) 


where dCij{s, .s)/d.§ is the parton luminosity for a channel involving two initial partons i,j with 
a center-of-mass energy y/§, and d'ij{s) is the partonic cross-section for that process. The sum 
is performed over all the partonic reactions that contribute to the background. An explicit 
expression for the parton luminosities in terms of the parton distribution functions is given in 


(A.l) in the appendix. 


Two simplifications can now be made. We assume the resonance to be narrow, hence the 


integral in (3.3) is performed over a small energy range As close to the resonance peak, with 
As <C The integrand can then be taken approximately constant over the integration range. 
If furthermore one performs the searches in a window of constant relative width As/m^, the 


equality (3.1) simplifies to 


2 f 2 \ - / 2 \ 2-^0 f 2 \'' / 2 ^ 

{iJ} 


{hi} 




(3.4) 


At energies much higher than the SM thresholds, the cross-sections dij decrease as s and 
one can finally write 


E' 

{*.i} 


dZj. 

ds 


L, 


\s,w?) = -j-Y^ 


dL 




L ds 

{*.i} 


{so,ml), 


(3.5) 


where the coefficients Cij = lim^-^oo saij{s) determine the composition of the background. While 
the s“^ behaviour of the cross-sections is clearly a good approximation for high-mass resonances, 
as was the case in [26| , one may wonder whether it can be used also in our context, where 
one is interested in scalar particles that are not too much heavier than the Higgs boson. We 


verified that (3.5) is actually valid down to a few hundreds of GeV, depending on the process. 


In any case, we do not use the experimental constraints below the region of validity for our 
extrapolations, and we do not provide extrapolated exclusions for masses below that limit. 
More details on this issue are given in appendix [A) 


Since the extrapolated values of m obtained from (3.5) are higher than the original mo for 


a collider with increased energy or luminosity, the range of masses where one gets a projection 
of the bounds will be shifted to higher values. Often, however, the physically more interesting 
region - at or below the TeV scale- lies outside that range, especially in the case of the colliders 
with highest luminosity or energy. In order to get an extrapolated exclusion in the full region 


of interest, in 26 the following method was used. One notices that, for a given energy of the 


new collider, the values of the masses where an exclusion is provided grows with the luminosity 
L. This follows immediately from (3.5) and from the fact that the parton luminosities are a 


decreasing function of s. Then, one can easily extend the extrapolation of the bounds to lower 
masses, decreasing the luminosity of the collider from L to L' < L; lowering the luminosity 
enough, one will eventually be able to get an exclusion for arbitrary masses down to the initial 
values. However, since the excluded cross-section increases with the inverse of L', the bounds 
for low masses will be much less powerful, and can be interpreted only as upper bounds on the 
actual reach of a collider. This is what has been done in . 

In order to get more realistic projections for low masses, in our analysis we rescale the cross- 
section limits, for each value of the mass, by the square root of the ratio between the nominal L 
of the new collider and the luminosity that was actually used for the extrapolation, V. This is 
a quite accurate approximation of the scaling of the excluded cross-section with the luminosity 


when the systematic errors are small (see (3.7)). 
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ATLAS 


CMS 


ZZ{U) 

y/s = 7-8 TeV, L = 4.6 + 20.7/fb 

39 


yT = 7 - 8 TeV, L = 5.1 + 19.7/fb 

40 

ZZ{2i2v) 

^/s = 7 TeV, L = 4.7/fb 

41 


55 

42 

ZZ{2i2j) 

5 ? 

43 


55 

42 

ZZ{212t) 

— 


55 

42 

WWi2£2i^) 

= 8 TeV, L = 21/fb 

44 


55 

42 

WW{£vjj) 

y^=7 TeV, L = 4.7/fb 

45 


y^ = 8 TeV, L = 19.3/fb 

46 

TT(all) 

— 



77 = 7 - 8 TeV, L = 5.1 + 19.7/fb 

42 

hh{4b) 

= 8 TeV, L = 19.5/fb 

47 


77 = 8 TeV, L = 17.9/fb 

48 

hh{blry^) 

Vs = 8 TeV, L = 20/fb 

49 


77 = 8 TeV, L = 19.7/fb 

50 


Table 2. Direct searches relevant for our analysis performed at the LHC experiments. 


To summarise, starting from an LHC exclusion [a x BR]o, the corresponding projection for 
a different collider, for a given mass m, reads 


\a X BRl (m: s, L) = min 


Let 


L'<L [y/LL' 


[a X BR]o(mo;so,Lo; 


mo(L') 


(3.6) 


where niQ^L') is determined from (3.5), substituting L with L'. 


Discussion of the main assumptions. Let us finally make a few comments on our prime 
assumption that the exclusion is driven by the background alone, which is not necessarily very 
accurate. First of all, the composition of the background must not change dramatically between 
the two experiments: the background that is dominant at the LHC must still be the leading one 
at higher energies. Next, the assumption is justified only if the analysis is done as a counting 
experiment of signal vs. background events, assuming the resonance to be narrow. If effects 
due to the shape of the resonance peak are taken into account, the results depend on signal 
and background kinematical distributions in a non-trivial wayj^ We always work in the narrow 
width approximation. The significance of a signal over a given background then is 


significance = 


Ns 


^NB + a^Nl + ^^Nl 


(3.7) 


where a and f3 are the systematical errors on the background and signal, respectively. In the 
hypothesis that the background dominates over the signal, which is always verified in our case, 
this last error can safely be neglected. However, in order to have the same significance at 
different collider experiments when Nb is hxed, one must assume that the systematic error a, 
when it is not negligible, remains approximately constant. 


A further issue is represented by the fact that the values of Ng and Nb in (3.7) in general 


depend on the acceptances and efficiencies, due to the experimental analysis and detector effects. 
They could vary for different experiments, especially when extrapolating over a wide range of 
e.g. from 8 TeV to 100 TeV. We will nevertheless assume, following 


26 , all the 


energies 

acceptances and efficiencies to be independent of the collider energy and the resonance mass. 
This assumption is another possible limitation of the method, but it is quite reasonable for 
the second stage of the LHC. For higher energy colliders, our results should be viewed as an 
estimate of what limits could be achieved assuming an experimental setup similar to the present 
one, and at the same time as a goal for the future detectors. 


^This is however what is done by the experimental collaborations at the LHC; we assume that the results do 
not deviate much from those obtained with a cut-and-count analysis. 
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Figure 3. Left: CMS 95% C.L. exclusion (black) on (Jpp^^ x from a combination of the ZZ 

and WW channels, together with the projections for the LHC at 13 TeV with 100 fb“^ (thin red), and 
at 14 TeV with 300 fb“^ (solid red) and 3000 fb“^ (red dashed). Right: reach in the hh —VV channel 
of the HL-LHC (dashed), HE-LHC (dot-dashed) and FCC-hh (dotted) with 3 ab“^. 

3.2 Current and future constraints 


Given that we are interested in models with a scalar singlet, the direct searches relevant to us 
are those for a heavy scalar decaying into W~^W~, ZZ and hh final states. The summary of all 
the searches that we have considered is provided in table These are only a partial list of all 
the searches for extra scalars carried out by LHC experiments, but they are the most sensitive 
ones for scenarios where the extra scalar is a singlet. 

Please note that all the extrapolations that we present here are valid for a generic spin-0 
resonance decaying into ZZ, WW and hh, and have therefore a general validity beyond the 
singlet scenario that we consider in this paper. 


The ZZ and WW channels. Many different searches for scalar resonances decaying into 
WW and ZZ m. fully leptonic and semileptonic final states have been presented both by ATLAS 
and CMS and are summarised in table The strongest individual constraint comes always from 
the ZZ —)• M and ZZ —)• 2l!.2v channels, which dominate respectively below and above about 
500 GeV. The WW and the other ZZ channels give a bound on the cross-section which is always 
weaker by a factor of a few. However, a recent analysis by CMS (4^ shows that the inclusion 


of all these final states is able to significantly improve the limit over all the mass rangej^ We 
will use only this global combination for our analysis in the di-vector channel. 

The background in the ZZ ^ M channel, which dominates the exclusion in the region of 
interest for our extrapolation, comes to a large extent from non-resonant ZZ production, and 
is mainly produced through the qq ^ ZZ partonic reaction [^. We therefore extrapolate 
this background using the Cqq parton luminosities, using everywhere the MSTW2008 parton 
distribution functions 


51 


The experimental exclusion is provided in the mass range from 145 GeV up to 1 TeV, but we 
consider only masses above 200 GeV because below that value the SM thresholds influence the 
shape of the background and the extrapolation procedure breaks down. We explicitly checked 
that above 200 GeV the excluded cross-section is driven by the parton luminosity to a very 
good approximation (see appendix [A]) . 

In figure]^ we show our results for the projected limits on x The left-hand 

panel shows the projections for the LHG at 13 TeV with an integrated luminosity of 100 fb“^, 
and at 14 TeV with a luminosity of 300 fb“^ and 3000 fb“^, together with the current constraint 
from GMS. The right-hand plot shows the extrapolations for the luminosity and energy upgrades 
of the LHG, respectively with a center-of-mass energy of 14 TeV and 33 TeV and both with a 
luminosity of 3 ab“^, and for the 100 TeV FGG-hh collider, also considering a luminosity of 3 


“In the same paper, the results are explicitly interpreted in the case of a scalar singlet extension of the SM. 
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Figure 4. Left: CMS 95% C.L. exclusion (black) on app^^ x BR^^hh from the 46 channel, together 
with the projections for the LHC at 13 TeV with 100 fb“^ (thin purple), and at 14 TeV with 300 fb“^ 
(solid purple) and 3000 fb“^ (purple dashed). Right: reach in the hh 46 channel of the HL-LHC 
(dashed), HE-LHC (dot-dashed) and FCC-hh (dotted) with 3 ab“^. 


ab“^. The reach of FCC-hh extends up to masses of about 11 TeV. 

We have checked that our limits are in very reasonable agreement with those obtained in 
the Snowmass study , where a full collider simulation of the signal and the background has 
been carried out, and with the ones for the HL-LHC presented by CMS in 


53 


One sees that the next phase of the LHC will improve the present bounds by more than a 
factor of 2 - with some improvement to be expected already in the first 13 TeV run - while 
another factor of 3 will be achieved by the high luminosity phase. In terms of the absolute limit 
on the cross-section, the high-energy colliders HE-LHC and FCC-hh will give weaker constraints, 
but one has to keep in mind that the cross-sections for a given mass rrifp will be much larger at 
such colliders. 


The hh channels. Searches for resonant hh production have been performed by both the 
ATLAS and CMS collaborations. The two most sensitive channels at present are hh —)• 46 
which gives the strongest exclusions for masses rrifp > 400 GeV, and hh —?■ 26 27 49,50 


47,48 


which is more constraining at low masses. Other decay channels, e.g. involving r leptons, could 


also become relevant for the next LHC run and future experiments 54 - 56 


The background in the 46 channel is dominated by tt and multijet events 48 . In terms of 


partonic processes, both gg and qq contribute, but the gluon-fusion background is expected to 
dominate at high energies, hence we will use the Cgg parton luminosities for the extrapolation. 
As in the previous case, we verihed that the inclusion of the other parton luminosities would 
change the result only marginally. We have checked explicitly that the extrapolation procedure 
is robust above about 500 GeV, where the ti threshold starts to influence the background (see 
appendix [a]), and therefore we restrict the exclusion to masses above this value. 

In figure we show our results for the projected limits on x from the 46 


channel, extrapolating the CMS exclusions [^. The energy and luminosity benchmarks and 
the notation are the same as in figure The improvements from LHC 8 to the future machines 
is roughly the same as for the VV channel. 

Concerning the 26 27 channel, its extrapolation is trustworthy only above 400 GeV, and in 
this region the 46 exclusion is always more powerful. For this reason we do not consider this 
channel for our future projections. We however include in our analysis the present exclusion 
which is relevant for low masses. We use also here the CMS data |50], since ATLAS sees a slight 
excess in this channel 


49 


3.3 Results and comparison of bounds 

We now interpret the bounds presented in the previous section in the context of a generic scalar 
singlet extension of the SM. 
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Figure 5. Reach of the direct searches for (j) — >• VV (red), (j) hh{4b) (purple), and (j) hh{2b2j) 
(magenta), for Vs = —75 GeV, and Xhs = As = 1. Grey: unphysical region. Left: current LHG 
exclusion at 95% G.L. (coloured regions), and projections for LHG13 (thin lines), LHG14 (solid lines) 
and HL-LHG (dashed lines). Right: comparison with the projections for HE-LHG (dot-dashed lines) and 
FCC-hh (dotted lines). 


The total production cross-section times branching ratio for (j) is given by (2.5) and (2.6). 


When dealing with current bounds we always use the SM Higgs total cross-section and branching 
ratios given by the Higgs Cross Section Working Group 57 , which are provided for masses up 


to 1 TeV. On the other hand, since the extrapolated limits extend up to much higher energies, 
where we cannot rely on the results of 57 , we use the LO quantities in our projections. In 
particular, we compute the gluon-fusion cross-section with the program higlu 58,59 , and we 
use MADGRAPH to calculate the vector boson fusion and VH associated production cross- 
sections. The sum of these three components constitutes a good approximation of the total 
Higgs production cross-section. Concerning the branching ratios, we use the results of for 
masses below 250 GeV, and we match them to the LO (plus NLO QCD) results above that 
scale, where the two quantities are practically equal. 

In figure a nd figure we show the present and future exclusions in the same niff, — Mhh 
plane of figure[^ for two different values of the singlet vev Vg- We fixed Xhs = A 5 = 1, however, 
as discussed above, the impact of their precise value on the results is very weak. The left-handed 
panels of both figures show the projections for the present and future stages of the LHC, up 
to the luminosity upgrade, while in the right-handed panels we show the range relevant for 
future high-energy colliders - HE-LHC and the 100 TeV FCC-hh. The constraints from the VV 
channel are drawn in red, and the ones from the two /i/i(46) and hh[2h2'y) channels in purple 
and magenta, respectively. As discussed above, we do not show future projections for the 2b2'y 
channel, since they are never relevant in the region where we perform the extrapolations. 

One sees immediately that, depending on the value of Vg-, the relative importance of the VV 
and hh searches can be very different. In particular, we show the cases of Vg = —75 GeV and 
Vg = 250 GeV as benchmarks of a large enhancement and suppression of BR^^hh, respectively. 
In the first case the searches in the hh channel dominate both at the 8 TeV and 14 TeV LHG, 
but for the larger masses that can be reached at future colliders their reach becomes comparable 
with the one from the VV channel. In the second case of figure on the other hand, 
is so strongly suppressed that searches in /i/i(46) will barely be sensitive to this scenario even 
at FGG-hh, while VV always gives a much stronger exclusion. 

The message that emerges from these figures is therefore pretty clear: at a high-energy 
collider, be it HF-LHG or FCC-hh, the VV searches will always put the strongest bound on 
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Figure 6. Reach of the direct searches in the VV (red) and hh{Ab) (purple) channels, for Vs = 250 GeV, 
and Xhs = = 1- The notation is the same as in figure 

the production of a singlet-like scalar. On the contrary, at the various phases of the LHC the 
hh —)• 46 process can be the most powerful probe, given some particular choice of the parameters. 
This is easily understood in terms of the fact that the enhancement (or suppression) of the hh 
branching ratio is particularly effective at low masses, where the deviations from the limit (2.7) 
can be larger (see (2.14)), but in the region of parameter space relevant for future hadron 
colliders (high masses and small mixing angle) it is not possible to modify more than 

a certain amount. Despite this strong dependence on the choice of parameters of the bounds in 
the VV and hh channels, the combined reach of the two exclusions is very similar in all cases: 
a suppression in will enhance the branching ratio into vectors, and vice versa, and the 

constraints on the cross-section in the two channels are of comparable order of magnitude. 

We further point out that the enhancement of is correlated with the modihcation 

of the triple Higgs coupling, as one can see from (2.14) and (2.15), and comparing figure]^ with 
figure This implies in particular that, if a resonance in the hh channel will be discovered in 
the future, in the singlet scenario under consideration one expects to see an enhancement also 
in non-resonant double Higgs production. 

In figure we show a comparison of all the direct searches with the indirect constraints 
from the measurement of Higgs signal strengths. Here we fix to the asymptotic value 

of 0.25, and the direct exclusion is always dominated by cj) ^ VV. As in figures [T] and 
we draw contours of constant instead of explicit indirect exclusions for each future collider, 
since the projections for the precisions that can be attained on the Higgs couplings are still quite 
uncertain. One can however easily read the numbers off tableNotice that in the right-handed 
plot the range of the vertical axis is changed in order to better show the region of small mixing 
angle most relevant for future colliders. 

The direct searches are always more powerful than the Higgs signal strength measurements 
for low values of the resonance mass, while the opposite is true for high masses. While this 
is a general feature of direct and indirect searches, such an interplay was a priori not guaran¬ 
teed in the physically interesting range of masses, and for the actual sensitivities of the future 
experiments. At the LHC, the point where the two searches become comparable varies from 

~ 500 GeV for LHC 8 to slightly more than 1 TeV at HL-LHC. Also a linear e'’'e“ collider, 
with a sensitivity to deviations in the Higgs couplings in the 1% ballpark, would become more 
sensitive than the HE-LHC upgrade above about 1 TeV (500 GeV for HL-LHC). The same mass 
of a TeV sets also the limit between FCC-ee and FCC-hh. 
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Figure 7. Left: Comparison between the combined direct LHC exclusion and reaches (red) and the 
deviations in Higgs couplings (pink), with = 0.25. Right: the same comparison, in the region 

relevant for future colliders. The coloured regions are excluded at 95% C.L. and the notation for the 
lines is the same as in figures and 

4 Non minimal supersymmetry 

We can now apply the previous general discussions to some concrete models that are particularly 
motivated scenarios for physics beyond the SM. In this section we consider the Next-to-Minimal 
Supersymmetric Standard Model (NMSSM), which is the minimal deformation of the MSSM 
that includes a scalar singlet in the spectrum. A chiral superfield 5, singlet under the SM gauge 
groups, is coupled to the two Higgs doublets Hu and of the MSSM through the superpotential 


Wnmssm = Wmssm + +/(5), (4.1) 

where f{S) is a generic polynomial up to third order in S. 

This deformation has important consequences for the naturalness of the theory, given the 
measured value of the Higgs mass. Indeed, in the MSSM a 125 GeV Higgs requires the largest 
contribution to nih to originate from SUSY-breaking effects, m\ < + A^, where is the 

radiative correction mainly due to the top-stop sector (see e.g. [61|). This in turn requires either 
top squarks in the multi-TeV range, or large trilinear terms, both cases implying a large tuning 
of the electroweak scale. The ultimate reason for this is that the supersymmetric couplings 
in the MSSM are controlled by the weak gauge couplings. The situation is different in the 
NMSSM, where a supersymmetric contribution to the Higgs mass originates at tree-level from 
the Yukawa-like term XSHuHd, 

ml<^^slp + mlclp + (4.2) 


where tan /3 = Vujvd is the ratio of the vev’s of the two Higgs doublets Hu^- As manifest from 
(4.2), the most natural scenario occurs when the coupling A is of order one, so that the largest 
contribution to rrih has a supersymmetric origin. We here briefly sketch the major features of 
such a picture: 


o It gives a large supersymmetric contribution to the Higgs mass (with a different depen¬ 
dence on tan/3 than the D-terms) that is sufficient to achieve 125 GeV at tree-level. 

o It modifies the dependence of the electro-weak vacuum expectation value on the La- 
grangian parameters, allowing for parametrically larger soft masses of stops and gluinos - 
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Figure 8. NMSSM with “strong” coupling, A = 1.2 and A = 70 GeV. Shaded regions and lines as in 
figure 

by a factor of 0{2\/g) - than in the MSSM, for a given amount of tuning [62f|64] . While 
a large A > 1 softens the tuning in u, it increases the tuning in the Higgs mass because it 
overshoots its measured value. 

o The coupling is not asymptotically free and if A > 0.7 the theory undergoes a strong 
coupling regime at a scale lower than the GUT scale [^. We will consider both regimes 
in the following, as they both have been deeply studied in the literature. 


The NMSSM with a sizeable A - which is the leading parameter also for phenomenology - opens 
then the interesting possibility that the lightest new particles be the extra scalar bosons of its 
extended Higgs sector (possibly with the exception of the LSP). These new physical degrees of 
freedom consist in three neutral scalars, two pseudoscalars, and one charged Higgs boson. 

The phenomenology of the NMSSM is controlled by several parameters, so that its study 
is often performed via scans on some motivated ranges (see e.g. [ 66 |), which need the spec¬ 
ification of a particular version of the NMSSM potential, i.e. the function f{S) and all the 
soft terms. Here we adopt a more general approach to have an analytic understanding of the 
phenomenological properties of the extra Higgs bosons, focussing on the CP-even ones, in the 


spirit of 28,67 


The SM-like Higgs held of (2.1) reads H = sin fi Hu + cos /3 Hd- As the present data suggest 


(see [M,67,^), a particularly interesting limit is the one where the orthogonal Higgs doublet, 
which does not participate in EWSB, is decoupled]^ Under the assumption of CP conservation, 
we can single out the CP-even component from the complex held S. In this case the system is 
equal to the two-state scenario described in section and can be studied along the same lines. 

The main parameters are again and with the former depending on the parameters 
of the model like ^ ^ 

^ n 2/3 -b A^. (4.3) 


Similarly to section]^ the mass eigenstates {h, (/>) are dehned as in (2.2), while the mixing angle 


The case of a non-decoupled second doublet is mostly probed via Higgs coupling measurements, and its 

Moreover, snch possibility is further constrained by 


phenomenology is similar to that of an MSSM Higgs 28 
the presence of a charged Higgs boson. 
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Figure 9. NMSSM with “perturbative” coupling, A = 0.7 and A = 80 GeV. Shaded regions and lines 
as in figure]^ 

7 between the Higgs and the singlet is now given by 


sin 7 = 


m\ 2/3 + (A^t;^/2) sin^ 2/3 + 




(4.4) 


The phenomenology is mainly set by 7 and m^f) but, as discussed in section [2T| quantities like 
the trilinear couplings depend on other parameters. For simplicity, in the following we will fix 
to its asymptotic value 1/4, so that there is no need to enter the details of the scalar 
potential!^ Notice however that all the analysis of section 2.2 applies, in general, also to the 
present case. In particular, from the point of view of phenomenology Vg is again the leading 
extra parameter. 


4.1 Results and discussion 

The quantity M‘^j^ now depends on parameters with a defined physical interpretation. Given 
the expression (4.4) we have chosen to fix (A, A) and let tan/3 free to vary. Notice that the 
value of A is almost irrelevant for sizeable A. We consider the following two cases: 


o “Strong” coupling, A = 1.2. This scenario is realised for example in A-SUSY 69 . The 


contribution from radiative corrections has been fixed to A = 70 GeV. This choice is 
consistent with a mass of the lightest stop both at the edge of the present exclusion (TOfTl 
for maximal mixing, and in the multi-TeV range. 


Yhe coefficients of the potential in (2.91 can be computed in terms of the parameters of the model. This 


amounts to have an explicit expression for f{S) in ( |4.1[ ) and for the soft term lagrangian Haoft. As an example 
in the scale invariant NMSSM (see for a review), 


W = XSHuHd+'^S\ 


-Aoft D AxiSHuHd + h.c.) + + h.c.) 


The matching for the coefficients, upon the decoupling of one of the two Higgs doublets, is given by 


A^ 


Xh='-j sm\2P) + g cos" (2/3) + 


A 5 ~ |/^| 5 A/fs = |A| , an = Axsm{2l3), as = 
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o Perturbative coupling, A = 0.7. In this case we use a value for A consistent with pertur¬ 
bative unification at the GUT scale, and we take A = 80 GeV in order to more easily 
reproduce the Higgs mass. 


Figures and show the present and future constraints on this scenario for both indirect 
and direct searches. We restrict ourselves to low values of tan/3 since this is the region most 
relevant for naturalness [62y63| . Gontrary to the generic Singlet case, the direct searches always 
dominate at the LHC and its upgrades. In the context of figure this can be interpreted 
with the observation that, for the physically motivated ranges of NMSSM parameters, M^h 
cannot be too large, see (4.3). The two benchmarks with different values of A have significantly 
different behaviours: the larger is A the smaller is the allowed parameter space. Moreover, while 
for a larger A leptonic colliders will start to be competitive with hadronic ones, for a smaller 
A only a per-mille level of precision in the Higgs couplings could be competitive with direct 
searches. This means that the observation of a deviation in the Higgs signal strengths, without 
any direct indication of an extra scalar, will strongly disfavour this scenario. It is interesting 
that, for supersymmetry at the LHC, our findings indicate that the strategies to best look for 
a singlet-like Higgs and a doublet-like one look quite different: direct searches in the first case, 
and mainly indirect ones in the second case (barring the large tan,0 case (7^[73] ). 

From the point of view of indirect searches, it is notable that also here large deviations from 
the SM value are possible in the triple Higgs coupling, as shown in 


28 and confirmed in a 


particular model by 74 , and contrary for example to the MSSM case. In general there are no 


particular differences with respect to the analysis of section given that Vg, A and Xs are free 
parameters also in the NMSSM. However, as we have already discussed, the NMSSM prefers 
low values of Mhh ■ Accordingly, the deviations from one in ghhh /dhM expected to be milder 

than in the generic case (see figures and [^ , but still allowing for observable effects at the 

LHC. 


5 Twin & Composite Higgs 


The presence of an additional scalar singlet is also typical of another class of models, namely 
those where the Higgs is a pseudo-Goldstone boson. In this type of models the Higgs boson 
emerges as an approximate Goldstone boson associated to the breaking of some global symmetry 
Q ^ In the case where there is only one physical pseudo Goldstone Higgs, the only other 
additional scalar singlet is the “radial mode”, i. e. the degree of freedom associated with vacuum 
expectation value /, with f > v, that breaks the global symmetry. Throughout this section we 
dub the extra physical scalar a instead of cj)- 

If the above description derives from a weakly coupled model at higher energies, then the 
radial mode a could be relatively light, significantly below 47r/. This is for example the case of 
the Twin Higgs (TH) idea [^, with a weakly coupled description 75- 
this possibility. 


79 , and we now turn to 


Twin Higgs. Twin Higgs models solve the little hierarchy problem up to some scale A, 
without requiring the presence of coloured degrees of freedom below such scale]^ This is contrary 
to the usual realisations of SUSY and composite Higgs models (CHM), where new coloured 
bosons or fermions are expected to lie within the LHC reach. In TH models the low energy 
degrees of freedom that cancel the “quadratic divergences” are total singlets of the SM and 
hence very difficult to detect at the LHC, while coloured particles are expected to be naturally 
heavier. In TH models even a null result from the LHC will leave the tuning of the weak scale 
at the level of 0(10%) (for recent studies see [79}|83] ). 

®Above A, a theory solving the big hierarchy between A and the Planck (or GUT) scale is needed. 
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In this paper we would like to emphasise that, if the TH has a weakly coupled description, an 
interesting possibility consists in trying to detect the “radial mode”, the mirror partner of the 
Higgs. In fact, differently from other mirror particles, it can be singly produced through a mixing 
with the Higgs. This is complementary to the idea that the only model-independent probe of 
such scenarios consists of the Higgs signal strengths. In order to discuss the phenomenology of 
this picture we here briefly review the basics of the model. 

o A copy of the SM gauge, matter and scalar content is added, related to the SM by a Z 2 
symmetry. A scalar portal coupling is allowed between the two sectors. 

o An important assumption is that the scalar potential is accidentally SO( 8 ) invariantin 
such a way we can have a spontaneous symmetry breaking SO(8)/SO(7) parametrized by 


= K {\Ha\^ + - y)' > 


(5.1) 


where Ha is the SM Higgs doublet, and 77 b is its Z 2 -symmetric partner, singlet under the 
SM gauge group. The global symmetry SO( 8 ) is manifest when the above expression is 


written in terms of the real components. From (5.1), and from the fact that the two copies 


of the SM are gauged, out of the 7 Goldstones only one remains as a physical fluctuation, 
the other six being “eaten” by the vector fields. In addition there is also a radial mode 
of mass ~ VX(/o. At this level one would like to explicitly break both SO( 8 ), to give a 
mass also to the other degree of freedom, and Z 2 , in order to get a hierarchy between the 
two vev’s in the two sectors. 

o This can be achieved by adding to the above potential other gauge-invariant terms, so 
that the full tree-level potential V reads 

V{Ha,Hb) = k{\Ha\^ + \Hb\^) + thWHaI^ - |77bP) + V,y^ra{HA, Hb) , (5.2) 

where the hrst term breaks SO( 8 ) but preserves Z 2 , and the second one breaks both 
groups. A nice feature of TH models is that the radiative corrections do not introduce 
any sensitivity to high scales, and the hnite contributions are suppressed by the fourth 
power of SM couplings (as an example k rsj ^(Ssm) 1^)13 


These three ingredients are sufficient to describe the phenomenology. The parameters in (5.2) 
are related to the vev’s (77^) = v‘^/2 and (77^) = (/^ — u ^)/2 via 


2 _ kA* Jq — {k + 2A*) 


(5.3) 


f2 ^ f2 

k(k + 2A*) ’ ■'°2A* + k’ 

so that / is the symmetry breaking scale of a strongly-coupled scenario, which indeed is obtained 
in the limit of large A*, i.e. when the mass of the radial mode is close to the cutoff and it can 
be consistently integrated out. 

From the diagonalisation of the mass matrix we can then fix nih and trade the residual free 
parameter for mg-, the mass eigenvalue of the singlet-like state. The mass eigenstates h and cr 


are related to the gauge ones via a rotation by an angle 7 as in ( 2 . 2 ), with 

P 


■ 2 

sm 7 = ^ 


mt 


mt - mt 


V 


(5.4) 


In the notation of section^ the parameter M?, is now given by 


Mhh = + 


(5.5) 


®For our discussion also SU(4)/SU(3) is fine. 

term k'\Ha\'^\Hb\'^ can be reabsorbed in a redefinition of k and of the coefficients in (5.11. The same is 


true for quadratic terms containing a different combination of [HAf and |i7s| 
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Figure 10. Model with BR^-^inv. = 3/7. 
(pink), excluded by direct searches (red), 
for the lines is as in figure 


Shaded regions: excluded at 95% C.L. by Higgs couplings 
• > rucr (blue), unphysical parameters (grey). The notation 


A crucial observation is that the only two free parameters are {mo-, /), which completely deter¬ 
mine the phenomenology. Like in the generic singlet case of section the mixing angle is also 
the value of the normalised couplings of h and a to SM particles. Other important informations 
are provided by the trilinear couplings of the mass eigenstates 


Shhh 
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2 r 


1 - 2v^/f 

y/l - uV/2 _ 


1 + 


3m,? / mf, 
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2m? 


Qcrhh — 
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2 r 
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1 - 


2m? (/^ — u^) 


+o(=4)l 

Vm?/ 


(5.6) 


Composite Higgs models. The above model only describes two real scalars, so that it 
closely resembles the case of a linearised minimal composite Higgs model with an SO(5)/SO(4) 
symmetry breaking [84] . In this context the global symmetry is spontaneously broken by some 
strong interaction, and the Higgs is again a pseudo-Goldsone boson of this breaking. The 
main phenomenological difference with the previous case is the presence of vector and fermion 
resonances related with the strong sector, with a mass of a few times the scale /□ 

There is one single radial mode S associated with this symmetry breaking pattern, which 
forms a quintuplet of SO(5) together with the Higgs doublet H] at energies higher than its 
mass one can describe this scenario through a linear cr-model. Although this might be an 
improper description of a strongly-interacting scenario, where the radial mode is generically 
expected to have a mass close to the cut-off of the theory dvr/, it provides a useful and simple 
parametrisation (38]. Using the same notation as above, and trading Hb with the pure real 


singlet S, the potential can be written as 86,68^ 


V{H, S) = A*(|Hp + 52 - flf + a/2|H|2 - (3\H?S\ 


(5.7) 


and it is not difficult to show that, to describe the scalar degrees of freedom, the potential in 


(5.7) is equivalent to the TH one in (5.2) (see e.g. footnote 10). 
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Figure 11. Model with BRo—>.inv. 
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= 0. Notation as in figure 


10 


5.1 Results and discussion 


Given the discussion above, the following phenomenological analysis will apply to the scalar 
sectors of both Twin (5.2) and composite (5.7) Higgs models. There are two main differences: 


o in the case of the “composite” Higgs, the fermion couplings of h (a) are not always a 
simple rescaling of the SM ones with cosy (sin 7 ), but they are model-dependent; 


o in a truly Twin Higgs model both h and a are expected to have some invisible decay to 
particles of the SM-copy. 

In the following we will assume a universal rescaling of the h and a couplings. It is however 
important to emphasise the differences between the two scenarios, commenting on the impact of 
an invisible decay channel for a in the TH modelj^ In the assumption that the two copies of the 
SM have the same gauge and Yukawa couplings, the expression for the two vev’s suggest that the 
particles of the invisible sector, which we indicate with a prime, have masses msM' — ^f/v. 

Given the approximate Z 2 symmetry, the decay widths of a to invisible particles (i.e. in the SM- 
copy) are similar to those of the SM Higgs boson to the corresponding visible particles. More 
precisely, the same predictions for the widths apply, upon substituting u 1 —>■ /. In particular, 
then, above the threshold ruo- > 2m\fi/i ~ 2xf/vxmw the invisible width can become dominant. 


ro--;.inv. — 3T„^Z'Z' — 3Tcr^zZ — ^^cr^hh, (5.8) 

where the last two equalities hold because of the equivalence theorem, that now involves the 7 


Goldstone bosons of our construction. This means that, instead of the relation (2.7), one has 

3 

7' 


BR„ 


3 HR 


(7—^ZZ 


(5.9) 


We will show the effect of the invisible branching ratios on direct searches, but we will also 
present results in the case BRq— i.inv. = 0, as a paradigmatic example of the composite Higgs 
picture. 


The results are shown in figures 10 and 11, where we present the case of a TH and a GH 
model respectively. The two scenarios differ mainly for the invisible branching ratio of the cr 


85 


and references therein. 


^^See e.g. 

^^In the case of TH, also h can have an invisible decay width, which is mainly due to decay to mirror bottoms. 


It is almost always under control for allowed values of /. We will not comment on this any further. 
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particle. Moreover, as already discussed, in the CH case we are assuming an universal rescaling 
of all the SM Higgs couplings, although this a model-dependent issue in concrete models. 

In both cases, if the constraints are presented in a plane (mo-, /) we see that at each collider 
stage indirect and direct searches will provide complementary information. In particular direct 
ones will help probing a region of small masses (and reasonably small mixing) that is not 
accessible to the Higgs coupling fit, allowing for another approach to study the TH models 
rather than the simple analysis of the deviation in the Higgs couplings. On the contrary the 
Higgs couplings put a robust bound for large values of m^-, where the direct searches are not 
really effective, and will be affected by the broadening of the resonance. Notice that the BRo— 


is always fixed by (5.6), differently from the generic and NMSSM cases, and implies that direct 


searches in TH and CHM are always dominated by the VV channel. 


5.2 Electroweak precision tests 


As already discussed 
T parameters 
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the IR-logarithms of (2.16) give sizeable contributions to the S and 
which depend on the mixing sin^ 7 and rua. In the case where the UV is 
mostly weakly coupled, hence for small values of the mass rria ~ /, EWPT are not an important 
constraint. 

On the contrary, in a strongly coupled scenario where m„ ~ dvr/, EWPT are expected 


to be more relevant. In this limit sin^y ~ while the logarithm in (2.16) is maximal, 


log(47r//m/j), and EWPT put a constraint on /. This reproduces the usual model independent 


contribution of composite Higgs models, and it can be read off (2.16). In the absence of other 
contributions, the EWPT bound becomes comparable to the Higgs coupling constraints, also 


given the signs of the S and T contributions of (2.16). In explicit composite models, tree level 


contributions from spin-1 resonances to the S parameter make the bound even stronger 87 
(see 


for a review). However, UV contributions can relax the bound from EWPT to a level 
comparable with the actual precision on Higgs couplings 89,90,26 , so we decide not to show 
this constraint, the result being more model-dependent. 


It is also interesting to compare the future improvements in the Higgs couplings and EWPT 
(see also 1^). EWPT are expected to benefit from a period of Z-pole measurements at an e"*'e“ 


collider, such as ECC-ee and CEPC. Present studies 91,92 suggest that the error on the S and 


T parameters (or the e-parameters) will shrink by an order of magnitude at a high-luminosity 
Z-pole facility. This improvement can be compared to the attainable precision on the hZZ 
coupling through the absolute measure of Uhz at a circular e^e“ collider, that is expected to 
strengthen the actual bounds by two order of magnitudes. They thus both reach a comparable 
level of precision on the scale / [26] : for this reason, and given the model dependency of the 
EWPT, we do not show projections of such constraints. 


6 Conclusions 


In this paper we have analysed the present and future constraints on models where a scalar 
singlet is added to the SM and is mixed with the Higgs boson. We focused on two main 
phenomenological aspects; direct searches for the new state at colliders, and indirect effects in 
the couplings of the 125 GeV Higgs due to its mixing with the singlet. 

We put a particular attention to the prospects of direct searches at the various stages of the 
LHC - the imminent 13-14 TeV run, and its possible luminosity and energy upgrades - and at 
a future high-energy hadron collider. In order to do that we have used an approximate method, 
already presented in p5, 26 , for extrapolating the existing constraints to higher energy and 


luminosity. The procedure is subject to several important assumptions and approximations: 
we have discussed them in detail and have performed various checks of our results, finding an 
overall agreement with the previous literature. 
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At the same time, we have considered the prospects for measuring deviations in the Higgs 
couplings, both at the LHC and at future linear and circular lepton colliders, trying to highlight 
the interplay between direct and indirect searches. 

We first discussed the model-independent case of a generic singlet, with interactions not 
constrained by any symmetry. The mixing of the SM Higgs boson with a scalar singlet allows 
an extremely simple description in terms of only two free parameters: the physical mass of the 
extra scalar and the mixing angle among the two states. These two quantities alone determine, 
in full generality, most of the phenomenology of the two particles, i.e. their production cross- 
sections and branching ratios into SM vectors and fermions. Interestingly enough, we have found 
that only one further parameter, the vev of the singlet Vg, is needed to describe the remaining 
properties of the Higgs system, namely the couplings among the scalars. Furthermore, the 
branching ratio of the singlet to a pair of Higgs bosons is fixed by the equivalence theorem at 
high masses, thus reducing effectively the total number of parameters to two in this limit. 

The relevant direct searches are always the ones for resonances in the W~^W~, ZZ and hh 
channels. Their relative importance depends on the size of BR^^hh, although the VV channel is 
typically the most sensitive. When compared with Higgs coupling measurements, direct searches 
will always provide a useful complementary tool; we find that, in general, they dominate for 
lower masses of the scalar singlet. An interesting point that emerged is the importance of the 
Higgs self-coupling: there are regions of the parameter space where large deviations from the 
SM could show up, even without observable effects in the Higgs signal strengths. This not only 
motivates an effort to optimise the reach at future machines, but also indicates that deviations 
could hide in the present and soon-to-come data, which we thus encourage to analyse in this 
respect. 

We have then specialised our analysis to a couple of interesting models that are expected 
to provide an extra scalar singlet at sufficiently low scales: the NMSSM and Twin Higgs and 
Composite Higgs models. These are at present also the best candidates for physics beyond the 
SM that are able to accommodate an (almost) natural electroweak scale with the various con¬ 
straints from the LHC. Their description deviates from the previous generic case since there are 
approximate or softly broken symmetries that constrain the interactions in the scalar potential. 
We then have found different limits with respect to the above case. 

1) The NMSSM with a light CP-even singlet. In this case direct searches provide an useful 
lamppost to find signs of new physics: we find it interesting that they are more powerful 
than the fit to the Higgs couplings. We considered the two cases of a perturbative and a 
strong coupling between the Higgs fields, the main difference among them is that the first 
case is more difficult to probe. Like in the case of a generic singlet, there is room for a 
sizeable modification of the trilinear Higgs coupling. 

2) Twin & Composite Higgs with a light “radial mode”. Although having a different physical 
interpretation, these two cases both provide a scenario where the Higgs boson is a pseudo- 
Goldstone boson of a spontaneously broken symmetry, and can be described in complete 
similarity. The real difference here is that the mixing angle does not go to zero with the 
mass of the singlet, providing a robust bound from the Higgs couplings on a large range 
of the parameter space. We have however found that if the mass of the radial mode is 
sufficiently small - a motivated assumption in a weakly coupled Twin Higgs scenario - 
then direct searches are stronger than the Higgs couplings. 

It is clear that the presence of an extra singlet scalar close to the weak scale is an interesting 
case to be studied. Despite its simplicity, it can offer diverse phenomenological scenarios, with 
interesting prospects for any future experiment both at lepton and hadron colliders. On the 
other hand, even without looking too far ahead, we have shown that already the second run of 
the LHC can efficiently probe this scenario, and we do hope that a great effort will be put by 
the experimental collaborations to close in on this picture. 
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A More on the extrapolation method 


In order to discuss the issues related to our extrapolation procedure, we would like to provide 
here more details about some of the steps already presented in section |3.1[ 


The estimation of the background relies on (3.3), where the main quantities are the SM 


cross-section aij and the parton luminosities d£ij/ds. The latter can be computed as 


93 


d£ 




ds 


(s,s) = 


1 + 




dx 


’ s/s 


X L 


fi{x,s)fj{ — ,s) + (i o j) 

\ ^ ^XS 


(A.l) 


where fi{x, Q'^) is the parton distribution function (PDE) of the parton i, and the factorisation 
scale has always been fixed to = s. Eor the background scaling with s to be driven by 
the parton luminosities only, one needs the cross-section to behave like 




(A.2) 


above a certain s value. Threshold effects due to SM particles decouple like m-gjyj/s, where msM 
is the mass of the particle relevant for the background, e.g. twice the top mass for pp —)• 46. 
When going above those thresholds, one could still be worried by the presence of logarithms, 
that are e.g. present whenever there is a contribution from a t-channel process. This is for 
example the case for qq —?■ ZZ, W~^W~, but also gg —?■ qq, which are all relevant backgrounds 
for the processes that we are interested in. Such logarithmic contributions are expected to 
be subdominant, because the experimental cuts used in the above searches tend to exclude 
the singular kinematic configurations which are responsible for them. Nonetheless, we have 
explicitly checked that the logarithms are not important at current energies, by verifying that the 
excluded cross-sections scale with the parton luminosities to an extremely good approximation 
(see figure 12) Einally, we verified that one can safely neglect the same logarithms also at 


higher energies, at least in the range relevant in this paper, by virtue of their weak dependence 
on s. 

We now discuss the parton luminosities used for the rescaling of the backgrounds relevant for 
our direct searches. The background for the process (p ^ ZZ, which leads the VV combination 


42 , is dominated at the partonic level by qq —)■ ZZ 40,42 . Matching (A.2) with the tree-level 
scattering cross-section qq ^ ZZ, one obtains that the coefficients for the up and down-type 
quarks satisfy the relation Cuujc^^ ~ ( 5 ^^ -k gnJ/iaid + dm) - where gLq and gRq are 
the Z couplings to the left- and right-handed quark currents. A further simplification can be 
achieved observing that, on the relevant mass scales at each collider, the parton luminosities of 
the up and down quark are in an approximate constant ratio. This means that the background 
qq —)■ ZZ can be rescaled to an extremely good approximation by considering either d£^j/ds or 


fact if the approximation of ( |A.2[ ) would not be verified, then one would see deviations in the above scaling 
due to the neglected logarithms. 
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Figure 12. Left: comparison between the LHC8 expected exclusion on app^cf, x BR^^vv (solid) and 
-\/L“^d£^j/ds (dashed). Right: the same as before, but with BR^^hh and £gg. 


dCuu/ds, given that in (3.5) any constant factor cancels out. Concerning now the other relevant 
search, the background to the process (j) —)• hh{4b) is consituted by tt and multijet events [is] , 
and pinning down its precise partonic composition is an involved task. However we expect an 
important contribution to come from the gg initial state, and we have chosen to perform the 
extrapolation using the gg parton luminosity. We have nevertheless checked that the impact of 
choosing another parton luminosity on the exclusions presented in this paper is mild. 

In figure [12] we show a comparison between the 95% C.L. expected exclusions on the cross- 
section times branching ratio, from LHC 8 , and the square-root of the (arbitrarily normalised) 
parton luminosities used for the corresponding extrapolations. The left panel shows the exclu¬ 
sion in the (p —?• VV channel together with the parton luminosity, while the right panel 
shows the exclusion in the p —)■ hh{4h) channel compared with the Cgg luminosity. One sees 
that, for masses above the relevant SM thresholds, the agreement between the two curves 
is extremely good. This confirms at the same time that the exclusion on the cross-section is 
determined mainly by the number of background events, and that the latter is proportional to 


the parton luminosities times a constant factor, thus substantiating equation (A.2). In other 


words, figure 12 provides a check that the approximation (A. 2 ) holds at the LHC 8 , and gives 


the further information of the mass from which it is valid. 
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